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We study the effect of electron-vibron interactions on the inelastic transport properties of single- 
molecule nanojunctions. We use the non-equilibrium Green's functions technique and a model 
Hamiltonian to calculate the effects of second-order diagrams (double-exchange DX and dressed- 
phonon DPH diagrams) on the electron-vibration interaction and consider their effects across the 
full range of parameter space. The DX diagram, corresponding to a vertex correction, introduces 
an effective dynamical renormalization of the electron-vibron coupling in both the purely inelastic 
and the inelastic-resonant features of the lETS. The purely inelastic features correspond to an 
applied bias around the energy of a vibron, while the inelastic-resonant features correspond to peaks 
(resonance) in the conductance. The DPH diagram affects only the inelastic resonant features. We 
also discuss the circumstances in which the second-order diagrams may be approximated in the 
study of more complex model systems. 

PACS numbers: PACS numbers: 71.38.-k, 73.40.-c, 85.65.-|-h, 73.63.-b 



I. INTRODUCTION 

Junctions consisting of a single organic molecule be- 
tween two metallic leads hold great promise for future 
nanoscale devices, where their potential applications in- 
clude switches, transistors, and sensors. Experimentally, 
it has proved difficult to control their production in an 
atomistic manner, and so theoretical studies are crucial 
for a full understanding of their behaviour. It is known 
that inelastic effects play an important role in the be- 
haviour of such deviceai"— , but as yet we lack a full un- 
derstanding of the processes at play that will lead to a 
complete interpretation of experimental results. 

In this paper, we use a model-system nanojunction in- 
cluding electron-vibration couplingii"— in order to de- 
termine, by investigating the whole parameter space, 
what level of diagrammatic expansion is appropriate to 
describe the electron-vibration interaction in such junc- 
tions. We cover the entire parameter space, thus account- 
ing for all the physical analogues to our model. The 
parameters we can vary correspond to the lead-molecule- 
lead coupling, the electron-vibron coupling strength, and 
the resonance of the electronic level with the leads' elec- 
tronic states 

In general, an organic molecule-based nanojunction is 
unlikely to have its HOMO or LUMO levels in align- 
ment with the equilibrium Fermi level of the leads, and 
so such a nanojunction will be dominated by what we 
term the ofF-resonant regime with strong tunneling at low 
bias. Moreover, coupling between the leads and the cen- 
tral molecule is likely to be relatively weak, in the sense 
that the corresponding coupling to the leads is much 
smaller than the corresponding hopping integrals in the 
leads themselves. Conversely, a system consisting of a 
nanoconstriction in a gold wire will have electronic levels 
in the constriction that are close to those of the leads, 
and a larger coupling between the central region and the 



leads, closer to the tight-binding hopping parameter of 
the leads. 

This work will distinguish between these regimes and 
discuss for which physical systems diagrams beyond 
the self-consistent Born approximation (SCBA) become 
relevant for electron-phonon coupling. We study this 
using a full non-equilibrium Green's-function (NEGF) 
tcchniqueii'i^'^'23i24^^j29-33j38^^ ^jjj(,j^ allows us to 

study all the different transport regimes in the presence of 
electron-vibron interaction. Following the spirit of many- 
body perturbation theory and Feynman diagrammatics, 
we include the commonly-used SCBA diagrams as well 
as second-order diagrams in terms of the electron-vibron 
interaction. A detailed description of the formalism and 
the numerical implementation of the NEGF code we have 
developed is given in Ref. jS^]. In this work we studied 
the equilibrium and non-equilibrium electronic structures 
of the nanojunctions in the presence of electron-vibron 
coupling. In the present paper, we now give and analyse 
results for the full non-equilibrium transport properties, 
namely the non-linear I{V) characteristics, the conduc- 
tance G{y) = dl/dV, and especially the lETS signal 
cPl/dV'^ calculated with our NEGF code. 

The paper is structured as follows. In Section [II] we 
summarize the key aspect of our methodology detailed in 
Ref. [53 . Results for the effects of the second order dia- 
grams on the non-equilibrium non-linear transport prop- 
erties are presented in Section IIIII They are separated 
into the features we observe for the purely inelastic effects 
in the lETS signal at bias equal to an integer multiple of 
the vibron energy, and for the inelastic resonant features 
related to the vibron replica of the electronic resonances. 
Our conclusions are given in Section [IVI In addition, we 
explain in detail in Appendix |X] how one of the second- 
order diagrams acts as a vertex correction to the Fock-like 
electron-vibron diagram. 



II. MODEL AND THEORY 

A fully atomistic description of the non-equilibrium in- 
elastic transport properties we wish to study is, unfor- 
tunately, beyond the reach of current ab initio methods. 
Instead we use a model system which retains the essential 
physics of the junction while reducing the calculations to 
a tractable size. A full description of both the model and 
our methodology is given in Ref. [57[, and so we review 
only the most salient features here. 

We use the single-site single-mode model (SSSM),in 
which the central molecule of the junction is represented 
by a single molecular level coupled to a single vibrational 
mode. We have already used this model and discussed 
its validity in our previous study on the equilibrium and 
non-equilibrium electronic structures of such a system 
coupled to two electron reservoirs^. 

The total Hamiltonian for the nanojunction is given by 



-fftotal — Hl + Hji + Vlc + VcR + Hr' + iJvib + H, 



^e-vib- 



(1) 



In this work, we represent the Hamiltonian of the left 
(L) and right (R) leads H^j^ with a non-interacting 
tight-binding model with semi-elliptic bands, although 
in principle it can take any valid form. The hopping be- 
tween leads and the central region is given by Vac = 
toa {Cad + d^Ca), whcrc toL,R IS the hopping integral be- 
tween the a = L,R lead and the central region. The 
central region contains the electron-vibron interactions. 
We choose that an electron couples linearly, via its den- 
sity, to the displacement of a single vibration mode. The 
Hamiltonian for the central region in the SSSM model is 
then given by 



Hn — Hr 



Hyih + He-vih 



~ Eod'd + hujQa^a + 7o(a + a)d^d, 



(2) 



where d^ (d) creates (annihilates) an electron in the 
molecular level Eq, which is coupled to the vibration mode 
of energy wq via the coupling constant 70 . 

A detailed analysis of the formalism of the non- 
equilibrium transport properties from NEGF and for in- 
teracting system is provided in Ref. [58|. The current 
la passing through each lead a is expressed in terms of 
two time Green's functions^. It is transformed into fre- 
quency representation for the steady-state regime to give: 

Ic = 11^ Tr{S]<H G>H - I]>H G<H}. (3) 

We vary the applied bias by moving the chemical po- 
tentials of the left and right leads. With the equilib- 
rium Fermi energies fi^'^ = ^°^ = /i°'^ = 0, we intro- 



duce a quantity rjv such that hl = M*"^ + rjveV and 
IJ,R = n"^ — {1 — riv)eV following Ref. [63. In this way 
we can create several forms for the potential drop across 
the junction. By setting rjy = 1, for example, we create 
an asymmetric drop whereby fiji remains constant while 




FIG. 1. The Fock-(F) and Hartree-like (H) diagrams 




FIG. 2. The dressed-phonon (DPH) and double exchange 
(DX) diagrams 



fiL is changed, whereas rjy ~ 0.5 gives a symmetric po- 
tential drop where /i^ rises (lowers) as /i/j lowers (rises) 
by the same amount. 

It now remains to construct our non-equilibrium 
Green's functions (details given in Ref. ISTI) . The retarded 
and advanced Green's functions are calculated using a 
Dyson equation 



C 



9c 



5^"E''^"G'^'°, 



(4) 



while the greater (G^) and lesser (G^) Green's functions 
arc obtained from a quantum kinetic equation with the 
form 

G>'< = (1 + G''S'-).g>^<(l + S^'G") + G''S>'<G^ (5) 

where gc is the non-interacting Green's function of the 
isolated central region. 

In this work we consider both first- and second-order 
contributions to the electron-vibron coupling. The first- 
order diagrams are shown in Fig.[Tl and, if calculated self- 
consistently, equate to the commonly-used self-consistent 
Born approximation (SCBA). Wc also make use of the 
two second-order diagrams, those which involve two 
phonon excitations (Fig. ^ . The first of these is similar 
in structure to the GW skeleton for electron-electron in- 
teractions, and consists of a Fock-like diagram where the 
phonon is dressed by a single electron-hole bubble (hence 
the appellation DPH for dressed-phonon diagram. Fig. ^ 
left). The second, which we call the double-exchange 
(DX) diagram (Fig. [Alright), includes two phonons simul- 
taneously, with the second being emitted before the first 
is reabsorbed. The DX diagram is part of the skeleton 
diagrams corresponding to vertex corrections. We use 
these diagrams to construct expressions for the electron- 
vibron self-energ3*S as the (total or partial) sum of each 
diagram: E^^^f' = Sf_,,, + S",,, + S^™ + 2°^,, . 

Note that in order to handle numerically sharply 
peaked functions or strongly discontinuous functions, we 



have found it necessary to include a very small but fi- 
nite imaginary part in our expression for the bare vibron 
Green's functioi*^. This also allows us to perform calcu- 
lations with a smaller number of w-grid points, as long 
as our imaginary part rj in the bare vibron Green's func- 
tion is around two to three times the w-grid spacing. We 
have already discussed in detail the effects of the corre- 
sponding extra broadening on the lincshapc of the spec- 
tral functions and on the values of the linear conductance 
in Ref. ^. 



III. RESULTS 

In this section we present the effects of the second- 
order diagrams on the full non-equilibrium transport 
properties of the nanojunction in the presence of electron- 
vibron coupling. Calculations of the Green's func- 
tions are performed with different levels of approx- 
imation for the electron-vibron self-energies (Figs. [T] 
and m. Fully self-consistent calculations using the 
first-order diagrams are annotated SCBA, those us- 
ing one or both second-order diagrams are annotated 
SC(BA-t-DX), SC(BA+DPH) or SC(BA+DX-KDPH) as 
appropriate. In addition we have performed non-self- 
consistent second-order corrections, i.e. by using the 
SCBA Green's functions to calculate the second-order 
diagrams, we then determine the new Green's functions 
without full self-consistency. These calculations are an- 
notated SCBA-hDX or SCBA-hDPH. 

The inelastic properties of the system are present in 
the current I{V) but are better represented by the second 
derivative of the current (Pl/dV^ as it is the signal that 
is directly measured experimentally in the form of the 
inelastic electron tunneling spectrum (IETS)i. 

The lETS curves present features, peaks or dips2^i^ 
at biases corresponding to the energy of a specific excita- 
tion, in our case to the energy of one or several excitations 
of the vibration mode. The peak feature is commonly 
associated with the opening of a new inelastic channel 
for the conductance of nanojunctions in the off-resonant 
regime, i.e. when the electronic level Sq is sufficiently far 
from the equilibrium Fermi level. In the case of the reso- 
nant transport regime (when eo is close to fi'^'^), a dip 
feature is obtained in the lETS. It is associated with 
electron-vibron backscattering effects and a decrease in 
the conductance at the threshold bias. Furthermore, be- 
ing the derivative of the conductance, the lETS curves 
also present features at biases corresponding to peaks in 
the conductance. We have found^ that in order to get 
a better aspect ratio for the lETS features correspond- 
ing to vibron excitations, it is more convenient to nor- 
malize the lETS curves by the dynamical conductance, 
i.e. [d'^I/dV^[dI/dV] = d/dV\nGiV) as in Refs.[l|,i- 
[1,0, i, and [13. 

We divide our results into two sections; the first for 
purely inelastic features, and the second for inelastic fea- 
tures associated with the electron resonance effects. The 
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FIG. 3. (Color online) The inelastic peak in the lETS around 
V — ujo = 0.3 for an asymmetric potential drop {-qv ~ 1.0, 
left panel) and a symmetric potential drop {rfv = 0.5, right 
panel). Results are shown for three different lead-molecule 
coupling parameters to ~ ioL.i? (black, red, and blue lines), 
and for two different electron-vibron couplings (medium 70 = 
0.195, solid lines and triangles; strong 70 = 0.25, dashed 
lines and circles). Lines represent calculations at the SCBA 
level, symbols SC(BA-(-DPH), from which we can see there 
is no distinguishable contribution to the inelastic peak from 
DPH for any parameter set. The other parameters are 
eo = +1-5, UJO = 0.3, 77 = 0.030. 



first category corresponds to features observed in the 
lETS signal at bias equal to an integer multiple of the 
vibron energy nujQ, i.e. a tunneling electron excites n vi- 
brons. The second category corresponds to inelastic reso- 
nant tunneling via the vibron replica associated with the 
main electronic resonance at Eq ^ eg — ^q/ujq, and hence 
are observed in the lETS signal for biases F ~ Eq ± "'^o- 
We will see below that the second-order diagrams have 
different effects on the lETS features depending upon 
their correspondence to one of these two categories. 



A. Purely inelastic features 

We consider first of all the off-resonant regime. In this 
limit the lETS features associated with inelastic resonant 
tunneling are sufficiently far (or sufficiently small for the 
higher vibron replica) from the inelastic feature dX V = 
Wo. We can thus avoid a superposition of the two different 
kinds of features. 

Results are presented in Fig. [3] for different electron- 
vibron couplings and for both symmetric and asymmetric 
potential drops. We note that, as would be expected intu- 
itively, the size of the feature increases with the electron- 
phonon coupling 70 — in fact the height with respect to 
the baseline is proportional to 79. As ioL,i? is increased, 
the lETS signal decreases in overall magnitude, although 
the feature at wq itself remains clear (Fig. [3] inset). 

In the off-resonant regime, we do not find any differ- 
ence to the curves when the second-order dressed-phonon 
(DPH) diagram is included. Possible reasons for this are 



discussed in section ITlI Bl 

We also note that the resuhs are dependent on the 
symmetry of the potential drops at the left and right con- 
tact. As shown in the separate panels of Fig.[3J the curves 
for symmetric potential drops have a much lower base- 
line than for the asymmetric potential drops, although 
the inelastic features have the same lineshape, position 
and magnitude (see Figs. [H |3] and [7]). This is because 
(for positive bias) in the asymmetric case, /i^ is rising 
bias while iir is kept fixed at its equilibrium value. Thus 
fj,L is approaching the electron resonance twice as fast as 
in the symmetric case, where fiL rises at the same rate 
as hr falls. The symmetric case therefore contains much 
less effect from the tail of the main electron resonance, 
and allows us to further isolate the purely inelastic part 
of the lETS signal. 

We now consider the effect of the double-exchange di- 
agram (Fig. [5] right) on the inelastic peak at ojq. Calcula- 
tion of S^^jjj is extremely computationally intensive, as 
it can be reduced neither to a simple convolution prod- 
uct, nor even to a simple double-convolution producl)^. 
Calculations for the integration in the energy representa- 
tion of the Green's functions and the self-energies scale as 
the cube of the number N^ of points in the energy-grid^. 
Moreover, as the finite imaginary part rj we have included 
in the vibron Green's function causes extra broadening, 
if taken too large it washes out much of the effect of 
the DX diagram. In order to decrease t] we need to 
increase the number of grid points, and thus, even for 
our minimal model, a fully self-consistent calculation of 
the double-exchange diagram becomes intractable with 
Ni^ > 2000. Hence for the DX calculations, wc usually 
work with N^ ~ 1500 - 2000, giving ?] - 0.03 - 0.04 for 
the total spectral width considered in our calculations. 

We work around this by calculating the effects of 
the DX diagram both in a self-consistent manner 
SC(BA-I-DX) or non-self-consistenly as a second-order 
correction to the SCBA result. 

The results of this are plotted in Fig.|4l We sec that the 
effect of the DX diagram is to both increase the height 
of the feature and to raise its baseline. For the set of 
parameters shown the self-consistency in the DX diagram 
calculations is not crucial. 

We see that the effect of DX can also be reproduced 
with a SCBA calculation in which we increase the value 
of 70- Therefore the DX diagram has the effect of renor- 
malizing 70 as it is part of the skeleton family of vertex 
correction, as shown in Fig. [5] We can thus approximate 
the effect of the DX diagram in the lETS with a Fock-like 
diagram with one renormalized vertex 70. 

The amplitude of the peak at ojo, instead of vary- 
ing as 7q, therefore now depends on 7070, and so we 
can define an effective clcctron-vibron coupling constant 
7o = VToTo j with 7o > 70- This allows us to make a 
more quantitative analysis of 70 by fitting the SCB A-t-DX 
curve to an SCBA curve with electron-vibron coupling 70, 
as shown in Fig. |4l 

Furthermore, we can study how 70 and 70 are cor- 
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FIG. 4. (Color online) Normalized lETS signal for the inelas- 
tic peak at ujq = 0.35. Results for SCBA and SC(BA-HDPH) 
are virtually identical. However, including the DX diagram, 
either self-consistently (solid red line) or non-self-consistently 
(red diamonds) raises the amplitude of the peak. This ef- 
fect can be approximated by an SCBA calculation done with 
a larger 70 (black triangles), indicating that DX generates 
an effective renormalization of 70. The other parameters are 
eo = +1.5, cjo = 0.35, ioL.H = 0.10, r? = 0.039, -qv = 0.5. 





FIG. 5. Second order double exchange (DX) diagram re- 
expression as a vertex correction to the Fock-like diagram, 
with a renormalized 70 — >■ 70 vertex. 



related by performing a series of calculations for differ- 
ent values of the parameters 70 and wq, and then fitting 
the SCBA-I-DX results to those from SCBA calculations. 
The results of this are shown in Fig. [SI from which we can 
see that, to a first approximation, the DX diagram con- 
sistently raises the effective electron-vibron interaction 
by approximately 3% within the range of parameters we 
used. 

However, we would like to point out that although the 
apparent effects of the DX diagram on the lETS sig- 
nal is to renormalize the coupling constant 70, the re- 
ality is much more subtle. In Appendix |^ we discuss 
in detail the renormalization effect of the DX diagram 
in terms of vertex corrections, and we show that such 
vertex corrections do not simply correspond to a map- 
ping of the SSSM Hamiltonian onto a similar Hamilto- 
nian with a static renormalization of the electron-vibron 
coupling constant, i.e. 70 — > 70 = 70 + |A|. Rather, the 
vertex correction actually generates a dynamical renor- 
malization of the electron-vibron coupling constant, i.e. 
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FIG. 6. (Color online) Effective static renormalized electron- 
vibron coupling parameter 70 as a function of the nominal 
coupling parameter 70 , from a fit of several sets of data com- 
paring SCBA+DX curves to SCBA (as in Fig. |4|. The data 
for all three values of a;o lie on a straight line with slope 1.03 
(the straight line with slope 1 is shown for comparison), imply- 
ing that the vertex correction increases the effective electron- 
vibron coupling by 3%. 



FIG. 7. (Color online) Normalized lETS signal for biases 
around bJa showing different SCBA fits for the renormal- 
ized electron-vibron coupling. The SCBA (dashed line) and 
SC(BA-(-DX) (solid red line) are shown for 70 = 0.34, together 
with SCBA calculations for 70 = 0.35 and 0.3525. Although 
these provide a good approximation, neither gives an exact fit 
to the lineshape of the SC(BA+DX) curve. The other param- 
eters are £0 = +1-5, 70 = 0.34, uo ~ 0.35,foL,fl ~ 0.10, r; = 
0.039, T]v = 1.0. 



70 ^7o(a;,a;')- 

This can be seen more clearly by considering SCBA 
calculations, for different values of the electron-vibron 
coupling parameter, and cheeking which of such calcula- 
tions correspond the best to a SC(BA+DX) calculation. 
The result is shown in Fig. jT] Although the difference 
between the best SCBA fits and the SC(BA-fDX) are 
not large, it is quite clear that a renormalized SCBA cal- 
culation does not provide exactly the same lineshape as 
a full SC(BA-I-DX) calculation for all the range of biases 
Y around the inelastic peak at wq- 

Finally we expect to see a peak feature at F ^ 2ujq 
in the lETS for the off-resonant transport regime. This 
peak feature is the two-vibration excitation equivalent of 
the feature observed at X^ = ojo- Since this a higher- 
order process, the amplitude of the feature should be 79 
times smaller than the feature at 1^ = wq. This feature 
has a rather small amplitude for all the electron-vibron 
coupling constants we have considered in this work since 
70 < 1. An example of a close-up of the lETS feature 
around V ^ 2ujq is given in Figure [8] We find that the 
amplitude of the peak with respect to the linear baseline 
is indeed approximately 79 = 0.112 (one order of mag- 
nitude) smaller than the corresponding amplitude of the 
peak at y = Wo (shown in Fig. |4]). Once more we find 
that the effects of the DPH diagram are negligible for 
this part of the lETS signal. In Fig. |S| we also include 
the results of a SCBA calculation for a larger coupling 
70 = 0.34 which mimics the renormalization effects of the 
DX diagram as discussion above. 
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FIG. 8. (Color online) Normalized lETS signal for biases 
around 2ajo for different levels of approximation: the SCBA 
(solid black line) and SC(BA-I-DPH) (red circles) are virtu- 
ally identical, as for the feature at wo- An SCBA calculation 
approximating the effect of DX is also included (dashed line). 
A straight baseline is included for clarity. The other param- 
eters are eo ~ +1-5, 70 ~ 0.335, (j;o = 0.35, ioL.fl ~ 0.10, r? = 
0.009, ryv =0.5. 



B. Inelastic resonant features 

The main electron-resonance peak occurs at the 
polaron-shifted value of eo ^ £o~7o/'^Oi ^-^d consists of a 
peak-dip feature in the lETS, as it corresponds to a reso- 



nant peak feature in the conductance. With no electron- 
phonon couphng (70 = 0), the lETS curve has no features 
other than the one corresponding to the resonant trans- 
mission in the conductance at T^ = eg- Once the electron- 
phonon couphng is turned on, phonon side-band peaks 
emerge in the spectral function at energies cj = Eq ± nujQ. 
The +nujQ features correspond to phonon emission (vi- 
bration excitation) by an electron, while the —nujQ fea- 
tures correspond to phonon-emission by a hole. In the 
lETS, they appear as peak-dip features (derivatives of an 
inelastic resonant peak in the conductance) with ampli- 
tude decreasing as the bias is further increased. At lower 
biases however, there are no features at w = Eq ~ 't-wq at 
the SCBA level except for very low values of the coupling 
to the leads {toL,R ^ 0.1, Fig. [9] top inset). 

Including the DPH diagram, however, introduces a 
small peak-dip features at Eq — wq just below the main 
peak-dip feature at V^ ^ eo in the lETS, as shown in 
Fig.inj). The DPH diagram also has a strong influence on 
the lineshape of the other phonon side-band peaks above 
Eo, increasingly so as Eq is brought within range of the 
equilibrium chemical potentials (i.e. eg ~^ in Fig. [S] 
bottom). 

We now consider the combined effects of the DPH di- 
agram as well as of the DX diagram on the specific case 
of the lETS feature at Eq — ujq. This is shown in Figure 
[TUl The DPH diagram increases strongly the peak-dip 
feature (at V ^ Eq — ujq) obtained from SCBA calcula- 
tions at medium/strong electron-vibron coupling. Note 
here the importance of the self-consistency in the calcu- 
lations: the second-order DPH diagram calculated as a 
second-order correction to SCBA (Fig. [TU] left) gives a 
completely wrong feature in the lETS. 

Interestingly, the self-consistent calculation with the 
DX diagram seems to give a similar feature to that ob- 
served in the SCBA calculations, but slightly shifted to- 
wards lower bias (Fig. [TUl right). This is completely con- 
sistent with the renormalization effects of the electron- 
vibron coupling by the DX diagram as discussed in the 
previous section. Indeed, the DX diagram renormalizes 
the coupling 70 towards a higher value 70. Consequently 
the renormalization of the molecular level by —j^/ujo is 
more important than for SCBA calculations, and thus 
the feature is moved towards lower bias. 

The calculations performed with both DX and DPH 
diagrams (Fig. [TU] right) generate a hybrid feature in the 
lETS in comparison to individual calculations with the 
second-order diagrams. However the new lETS is not 
simply obtained by a linear superposition of the individ- 
ual effects of the DPH and DX diagrams. 

It might at first seem strange that the DPH self-energy 
is ncgligblc at V = nujQ, where one might expect it to 
be influential, but that it has a significant effect at bi- 
ases > £o- While this is to some extent related to the 
strength of the electron-vibron coupling, there is another, 
more important, underlying cause. The DPH diagram in- 
volves an electron-hole bubble, and so for this diagram 
to become relevant, simultaneous electron and hole states 
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FIG. 9. (Color online) Normalized lETS signal. Top: In- 
elastic resonant features for different values of 70. With no 
electron-vibron coupling (solid black line), there is a single 
feature at eo = 0.7. With electron-vibron coupling (red and 
blue dashed lines), the main peak moves to £0 ~ £0 — 7o/'^o 
and side-band peaks appear. The inset shows a tiny feature 
at eo — a;o that only occurs for to ^ 0.1 in SCBA calcula- 
tions but reappears for all to once DPH is included. Bottom: 
The inelastic resonant features for various eo at 70 = 0.15 
for both SCBA (lines) and SC(BA+DPH) (symbols) calcu- 
lations, showing the increasing influence of DPH in the res- 
onant transport regime eo — i- 0. The other parameters are: 
to = 0.2, ujo = 0.3, riv = 1.0, 77 = 0.03. 



must be available. This is not the case when the spec- 
tral functions of the coupled electron-vibron system are 
mostly empty or mostly filled. When the bias is sig- 
nificantly low and both Fermi levels ^j-l.r are below the 
electron resonance level, these excitations are inaccessible 
and so there is no effects from the DPH diagram. Once 
the bias is increased to within range of Eq, or £0 i nwo 
however, these electron-hole states become accessible and 
the DPH diagram becomes influential (unless the spec- 
tral function is mostly filled). This is also borne out by 



>7 

o 

— 

CO 

^ ft 



5- 



1 ' 1 ' 


' 1 




9 


SCBA 




-- SCBA+DPH 








■ -■ SC(BA+DPH) 






8 






f 






1/ 




- 




l/' - 


7 


^■''/""v J 


V 




1 /( 


/' 






1 A\ 


I 


- 


6 


■ // 








V/ \ 


1 






# 


1 






yy 


\ I 




3 


1 1 1 1 


, 1 







I ' r 

SCBA 

SC(BA+DX) 

SC(BA+DX+DPH) 




0.75 0.8 



0.85 0.9 
V 



0.95 



0.75 



0.8 0.85 0.9 
V 



0.95 



FIG. 10. (Color online) Normalised lETS signal for biases 
around the first vibration side-band peak 1/ ~ eo — ojo = 
1.19 — 0.4 ~ 0.8. Calculations performed for different ap- 
proximations for the electron-vibron self-energy: (left panel) 
SCBA (black solid line), SCBA plus second order DPH correc- 
tion (dashed red line), SC(BA+DPH) (dot-dashed blue line); 
(right panel): SCBA (solid black hue), SC(BA-fDX) (dashed 
red line) and SC(BA+DX-^DPH) (dot-dashed blue line). The 
other parameters are eo = +1-5, 70 = 0.35, a;o = 0.4, ioL.H = 
0.09, 7? = 0.039, T?v = 1. 



the increasing contribution from DPH to the hneshape of 
the phonon side-band peaks as the electron level eois de- 
creased, moving DPH's sphere of influence to lower and 
lower biases (see lower panel in Fig.[S]). 



C. Summary over the entire parameter range 

Having examined the role of the second-order diagrams 
in detail for characteristic selected sets of parameters, we 
now show results across the entire parameter range. In 
order to present this in a concise manner, we have com- 
piled maps of our lETS results comparing SCBA calcu- 
lations to those with SC(BA-I-DPH) indicating in which 
regimes the DPH self-energy has a significant effect. For 
the off-resonant regime (Fig. [TT] left), we can see that the 
greatest effect of DPH is apparent at higher bias (i.e. ap- 
proaching the electron resonance Eq) and when the lead- 
molecule-lead coupling is small. For the resonant case, 
however, we can sec that the DPH self-energy only gives 
a non-negligble contribution in the region of parameter 
space where the coupling to the leads is small, and at low 
bias. 

This has potential implications for real molecular junc- 
tions. Consider a junction which has its dominant molec- 
ular levels far from the equilibrium Fermi levels of the 
leads. At sufficiently high bias, the DPH self-energy will 
have a significant contribution to the inelastic spectra. 
This could occur in the case of a junction formed from 
an organic molecule if the electronic level of the molecule 
is within range of the intended operational bias of the 
junction. If however the dominant molecular electronic 
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FIG. 11. (Color online) Maps generated from multiple cal- 
culations of SCBA and SC(DA+DPH) lET spectra for the 
off-resonant (left panel) and resonant (right panel) regimes 
showing the normalized absolute ratio of the SCBA and 



SC(BA-hDPH) spectra 1 



-'SCCBA + DPH) 



It gives a near- 



lETSsCBA 

zero result (blue and green areas) when the two calculations 
give the same spectrum, and a positive number (red areas, 
i.e. bottom right corner for SCBA and bottom left corner 
for SC(BA-I-DPH)) where there is a substantial difference be- 
tween the two. The other parameters are eo = 1.5/0.0 (off- 
resonant/resonant), 70 = 0.195,0^0 = 0.3. 



level is close to the leads' Fermi levels, and the coupling 
to the leads is large (as would be the case, for exam- 
ple, in a gold nanoconstriction) then DPH will not give 
a significant contribution at any applied bias. 



IV. CONCLUSIONS 

By using the non-equilibrium Green's functions tech- 
nique, we have studied the effect of electron-vibron in- 
teraction on the inelastic transport properties of single- 
molecule nanojunctions for a model system. We have 
included not only the first-order diagrams (BA) but 
also the second-order diagrams (double-exchange DX 
and dressed phonon DPH diagrams) for the electron- 
vibration interaction. We have calculated the inelastic 
electron tunneling spectrum (lETS) across the full range 
of parameters available to our model. The effects of the 
second-order DX and DPH diagrams are different and 
affect different features of the lETS signal. The effects 
of these diagrams are generally less visible in the inte- 
grated quantities, such as the current or the derivated 
lETS signal, than in the spectral functions^. However 
their effects are non- negligible, and are important for the 
full understanding of the spectroscopic information con- 
veyed by the lETS signal. 

The effect of the dressed-phonon (DPH) diagram is 
more important in the bias regions where one of the leads' 
chemical potentials begins to impinge upon the electron 
resonance or one of its vibron replica (i.e. for resonant 
inelastic features). Its effect is reduced both by increas- 
ing the lead-molecule-lead coupling and/or reducing the 
electron-vibron interaction. The renormalization of the 
vibron propagator (DPH) has been shown to be strongly 
dependent on the self-consistency of the calculations. It 
would be interesting now to study the effects of the full 
series of the electron-hole bubble on the renormalised vi- 



bron propagator (i.e. full GVF-likc diagram). 

The double-exchange diagram (DX) affects all the fea- 
tures in the lETS signal (i.e. resonant inelastic features 
and purely inelastic features at V '^ ujq)- The correc- 
tions are small in the weak-to-medium electron-vibron 
coupling because they are of the order of 0(7g). How- 
ever we have shown, numerically and analytically, that 
the effect of DX is similar to a dynamical renormalization 
of one vertex in the Fock-like diagram. More interest- 
ingly, the complex form of the non-equilibrium dynamical 
renormalizcd electron-vibron coupling 70(0;, w') we have 
derived analytically can be adequately replaced in our 
lETS calculations by a single static renormalized param- 
eter 7o. This important result leads us to believe that 
the second-order DX calculations, which are extremely 
costly in computing time even for our model system, can 
be incorporated in calculations for realistic systems by 
an appropriate renormalization of the vertex in a low- 
computing-cost SCBA calculation. 



Appendix A: Non-equilibrium vertex corrections to 
the Fock diagram 

In Appendix A of Ref . [53| we have given all the details 
for the derivations of the first- and second-order electron- 
vibron self-energies. In this section we show how the 
second-order double-exchange (DX) diagram can be re- 
cast in an effective first-order Fock-likc diagram with a 
renormalized vertex. 

We recall that for the SSSM model, the Fock and DX 
self-energies defined on the Keldysh contour Ck are given 

by 



The above expression for the vertex function is com- 
patible with the second-order expansion of the electron- 
vibron interaction. A generalization of the vertex func- 
tion (see Fock-like diagram in Figure [5]) to all orders 
of the interaction is possible, though beyond the scope 
of the present paper. Note that at the lowest or- 
der, the renormalized electron-vibron coupling parameter 
would be given by 7o(T3,r4;T2) = 7or''°'('r3, 7"4; ^2) with 
r(0)(r3,T4;T2) = 6{t3 - T2)6{t^ - T2). Hence Eq.(|X2]) 
would simply be transformed Eg. ljAip as expected. 

Using the rules of analytical continuation on the real- 
time branches given in Appendix A of paper Ref. |57| , we 
find the different components of the self-energies. Then 
after taking the Fourier transform of the different quanti- 
ties in the steady-state limit, i.e. X{t, t') = X{t — t'), we 
find the following expression for the different components 
of the energy-dependent self-energies: 



Sfii^^^^ (^) = n'o I ^ G'^''^' («) D^o"^' i^-v): (A6) 



and 



C3,C4 



r,C3C4C2 



l'o'''''{u;,v), 



(A7) 



where the non-equilibrium dynamical renormalised 
electron-vibron coupling is given by 



Sf.vib(n,r2) == i7o' Doir,,r2) G(Ti,r2), (Al) ^^^'^^^^co , v) ^ 'jl f ^ G'^^i^iv - u) 

3 J 27r 



and 



(A8) 



Sfv'[b(n,T2) = -^ / dT3dT4G(Tl,T3)I?o(Tl,T4)x 
^ JCk 

G(t3,T4)Do(t3,T2)G(t4,T2). 

(A2) 

The DX self-energy can be rewritten as an effective Fock- 
like diagram after introducing a renormalized electron- 
vibron coupling parameter 7o(T3,r4;r2) and the vertex 
function F (r3 , T4 ; T2 ) : 



D^„''^'{u)G'^-'^''{uj-u). 



Sf ^b (n , T2 ) = i7o / dr3 dr4 70 (T3 , r4 ; r2 ) x 

JCk 



Do{ti,T4)G{ti,T3), 



with 



i^ 

loirs, T4;t2) = -^r(r3, T4; T2), 



(A3) 



(A4) 



and where the vertex function is given by 

r(T3,T4;r2) = G(T3,r4)i:)o(T-3,T2)G(T4,T2). (A5) 



(The index Ci = ± labels the branch of the Keldysh 
time-loop contour Ck and are related to the usual con- 
vention: time-ordered (t ~ -\ — h), anti time-ordered (t = 
), greater (>= — +) and lesser (<= H — ) compo- 
nents.) 

Because all the quantities are originally defined on C'k 
and because the vertex function is a 3-point (3 times) 
function, the non-equilibrium dynamical renormalized 
electron-vibron coupling is a complex function of three 
indices Q and of two energy variables. Such a dynamical 
renormalization (including non-equilibrium conditions) is 
much more complicated than a simple static renormaliza- 
tion of the electron-vibron constant coupling 70 ^^ 7o- 

In our analysis of the lETS signal in the off-resonant 
transport regime, we try to keep the interpretation of 
the results as simple as possible, and we show that the 
renormalization of the lETS signal due to the DX dia- 
gram can be fairly well approximated by a simple static 
renormalization of the coupling constant 70 for applied 
bias around the vibration frequency V ^ ujq ± 20%. 



ACKNOWLEDGMENTS 



This work was funded in part by the European Com- 
munity's Seventh Framework Programme (FP7/2007- 



2013) under grant agreement no 211956 (ETSF e-I3 
grant). 



* |louise.dash@york.ac.uk| 

^ K. W. Hipps and U. Mazur, Journal of Physical Chemistry 

97, 7803 (1993). 
^ N. Liu, N. A. Pradhan, and W. Ho, Journal of Chemical 

Physics 120, 11371 (2004). 
^ J. G. Kushmerick, J. Lazorcik, C. H. Patterson, and 

R. Shashidhar, Nano Letters 4, 639 (2004). 

* L. H. Yu, Z. K. Keane, J. W. Ciszek, L. Cheng, M. P. 
Stewart, J. M. Tour, and D. Natelson, Physical Review 
Letters 93, 266802 (2004). 

^ L. H. Yu, C. D. Zangmeister, and J. G. Kushmerick, Nano 

Letters 6, 2515 (2006). 
« D.-H. Chae, J. F. Berry S. Jung, F. A. C. C. A. Murillo, 

and Z. Yao, Nano Letters 6, 165 (2006). 
'' J. M. Beebe, H. J. Moore, T. R. Lee, and J. G. Kushmerick, 

Nano Letters 7, 1364 (2007). 
^ N. Okabayashi, Y. Konda, and T. Komeda, Physical Re- 
view Letters 100, 217801 (2008). 
® H. Gawronski, J. Fransson, and K. Morgenstern, Imaging 

of inelastic waves in iets maps (2009), arXiv:0911.4053vl. 
^° Y. Kim, H. Song, F. Strigl, H.-F. Pernau, and T. L. and. 

Mechanical control of vibrational states in single-molecule 

junctions, arXiv:1011.3226vl. 
" P. Hyldgaard, S. Hershfield, J. H. Davies, and J. W. 

Wilkins, Annals of Physics 236, 1 (1994). 
^^ H. Ness and A. J. Fisher, Physical Review Letters 83, 452 

(1999). 
^^ H. Ness, S. A. Shevlin, and A. J. Fisher, Physical Review 

B 63, 125422 (2001). 
^"^ H. Ness and A. J. Fisher, Europhysics Letters 57, 885 

(2002). 
^^ K. Flensberg, Physical Review B 68, 205323 (2003). 
" T. Mil, S. Tikhodeev, and H. Ueba, Physical Review B 68, 

205406 (2003). 
^^ M. J. Montgomery, J. Hoekstra, A. P. Sutton, and T. N. 

Todorov, Journal of Physics: Condensed Matter 15, 731 

(2003). 
^® A. Troisi, M. A. Ratner, and A. Nitzan, Journal of Chem- 
ical Physics 118, 6072 (2003). 
" Y. C. Chen, M. Zwolak, and M. di Ventra, Nano Letters 

4, 1709 (2005). 
■^° N. Lorente and M. Persson, Physical Review Letters 85, 

2997 (2000). 
^^ T. Frederiksen, M. Brandbyge, N. Lorente, and A. P. 

Jauho, Physical Review Letters 93, 256601 (2004). 
^^ M. Galperin, M. A. Ratner, and A. Nitzan, Nano Letters 

4, 1605 (2004). 
^^ M. Galperin, M. A. Ratner, and A. Nitzan, Journal of 

Chemical Physics 121, 11965 (2004). 
^* A. Mitra, I. Aleiner, and A. J. MiUis, Physical Review B 

69, 245302 (2004). 
^^ A. Pecchia, A. di Carlo, A. Gagliardi, S. Sauna, T. Frauen- 

hein, and R. Gutierrez, Nano Letters 4, 2109 (2004). 
^^ A. Pecchia and A. di Carlo, Reports on Progress in Physics 

67, 1497 (2004). 



^^ Z. Chen, R. Lii, and B. Zhu, Physical Review B 71, 165324 
(2005). 

^^ M. Paulsson, T. Frederiksen, and M. Brandbyge, Physical 
Review B 72, 201101 (2005). 

^^ D. A. Ryndyk and J. Keller, Physical Review B 71, 073305 
(2005). 

^^ N. Sergueev, D. Roubtsov, and H. Guo, Physical Review 
Letters 95, 146803 (2005). 

^^ J. K. Viljas, J. C. Cuevas, F. Pauly, and M. Hafner, Phys- 
ical Review B 72, 245415 (2005). 

^^ T. Yamamoto, K. Watanabe, and S. Watanabe, Physical 
Review Letters 95, 065501 (2005). 

^^ A. Cresti, G. Grosso, and G. P. Parravicini, Journal of 
Physics: Condensed Matter 18, 10059 (2006). 

^^ M. Kula, J. Jiang, and Y. Luo, Nano Letters 6, 1693 
(2006). 

^^ M. Paulsson, T. Frederiksen, and M. Brandbyge, Nano Let- 
ters 6, 258 (2006). 

'■'^ D. A. Ryndyk, M. Hartung, and G. Cuniberti, Physical 
Review B 73, 045420 (2006). 

^^ A. Troisi and M. A. Ratner, Nano Letters 6, 1784 (2006). 

^* L. de la Vega, A. Martm-Rodero, N. Agra'it, and A. Levy- 
Yeyati, Physical Review B 73, 075428 (2006). 

^® M. C. Toroker and U. Peskin, Journal of Chemical Physics 
127, 154706 (2007). 

■"^ T. Frederiksen, M. Paulsson, M. Brandbyge, and A.-P. 
Jauho, Physical Review B 75, 205413 (2007). 

"^ M. Galperin, A. Nitzan, and M. A. Ratner, Physical Re- 
view B 74, 075326 (2007). 

"'^ D. A. Ryndyk and G. Cuniberti, Physical Review B 76, 
155430 (2007). 

■*^ B. B. Schmidt, M. H. Hettler, and G. Schon, Physical Re- 
view B 75, 115125 (2007). 

'''' A. Troisi, J. M. Beebe, L. B. Picraux, R. D. van Zee, D. R. 
Stewart, M. A. Ratner, and J. G. Kushmerick, Proceed- 
ings of the National Academy of Sciences of the USA 104, 
14255 (2007). 

"^ Y. Asai, Physical Review B 78, 045434 (2008). 

**^ C. Benesch, M. Ci'zek, J. Klimes, L Kondov, M. Thoss, and 
W. Domcke, Journal of Physical Chemistry C 112, 9880 
(2008). 

'''^ M. Paulsson, T. Frederiksen, H. Ueba, N. Lorente, and 
M. Brandbyge, Physical Review Letters 100, 226604 
(2008). 

■'^ R. Egger and A. O. Gogolin, Physical Review B 77, 113405 
(2008). 

"^ S. Monturet and N. Lorente, Physical Review B 78, 035445 
(2008). 

^° E. J. McEniry, T. Frederiksen, T. N. Todorov, D. Dundas, 
and A. P. Horsfield, Physical Review B 78, 035446 (2008). 

""^ D. A. Ryndyk, R. Gutirrez, B. Song, and G. Cuniberti, in 
Energy Transfer Dynamics in Biomaterial Systems, edited 
by I. Burghardt, V. May, D. A. Micha, and E. R. Bittner 
(Springer Berlin Heidelberg, 2009), vol. 93 of Springer Se- 
ries m Chemical Physics, pp. 213-335, ISBN 978-3-642- 



10 



02306-4. 
^^ B. B. Schmidt, M. H. Hettler, and G. Schon, Physical 

Review B (Condensed Matter and Materials Physics) 77, 

165337 (2008). 
^^ M. Tsukada and K. Mitsutake, Journal of the Physical 

Society of Japan 78, 084701 (2009). 
^* J. Loos, T. Koch, A. Alvermann, A. R. Bishop, and 

H. Fehske, Journal of Physics: Condensed Matter 21, 

395601 (2009). 
'"^ R. Avriller and A. L. Yeyati, Physical Review B 80, 041309 

(2009). 
^•^ F. Haupt, T. Novotny, and W. Belzig, Phys. Rev. Lett. 

103, 136601 (2009). 
" L. K. Dash, H. Ness, and R. W. Godby, Journal of Chem- 
ical Physics 132, 104113 (2010). 
^* H. Ness, L. Dash, and R. W. Godby, Physical Review B 



82, 085426 (2010). 
^'■^ H. Wang and M. Thoss, Numerically exact, time- dependent 

treatment of vibrationally coupled electron transport in 

smgle-niolecule junctions (2011), arXiv:1103.4945vl. 
^° A. Garcia-Lekue, D. Sanchez- Portal, A. Arnau, and 

T. Frederiksen, Simulation of inelastic electron tunneling 

spectroscopy of single molecules with functionalized tips 

(2011), arXiv:1103.4302vl. 
^^ A. Ueda, O. Entin-Wohlman, and A. Aharony, Effects of 

coupling to vibrational modes on the ac conductance of 

molecular junctions (2011), arXiv:1101.4440vl. 
•^^ Y. Meir and N. S. Wingreen, Physical Review Letters 68, 

2512 (1992). 
®^ S. Datta, W. D. Tian, S. H. Hong, R. Reifenberger, J. L 

Henderson, and C. P. Kubiak, Physical Review Letters 79, 

2530 (1997). 



